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^ ■ Abstract 

o 

T^lj- , Quantum circuit network is a set of circuits that implements a certain computation task. Being 



o 



at the center of the quantum circuit network, the multi-qubit controlled phase shift is one of the 
most important quantum gates. In this paper, we apply the method of modular structuring in 
classical computer architecture to quantum computer and give a recursive realization of the multi- 



Q^, qubit phase gate. This realization of the controlled phase shift gate is convenient in realizing certain 



quantum algorithms. We have experimentally implemented this modularized multi-qubit controlled 

X 

;h ' phase gate in a three qubit nuclear magnetic resonance (NMR) quantum system. The network is 

demonstrated experimentally using line selective pulses in NMR technique. The procedure has the 
advantage of being simple and easy to implement. 
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I. INTRODUCTION 

Quantum computer is born with the combination of quantum mechanics and computer 
science. Quantum properties have enabled quantum computer to outperform a classical 
computer for factorizing a large number [ll , or searching an unordered database [j] • In Ref . 
[3|, Deutsch proposed quantum gates and quantum network, and described a blue print for 
quantum computer. It has been shown that it is possible to construct arbitrary n-qubit 
quantum gate by using only a finite set of one-qubit gates and two-qubit gates |^, |5| . These 
basic quantum gates are universal for quantum computation[6j. Barenco et al, and Cleve 







have developed methods for designing networks in multiple- qubit system 

Quantum network can be viewed as the quantum analog of the classical computing. It 
provides with us a convenient method to build a quantum computer similar to building a 
classical computer. A classical compute can be implemented as a fixed classical gate array, 
with an input programm and data. A universal gate array can be programmed to perform 
any possible function on the input data. However for a quantum computer, Nielsen and 
Chuang demonstrated a rather different property [9]. They pointed out it is not possible 
to build a fixed, general purpose quantum computer which can be programmed to perform 
any arbitrary quantum computation task. The implementation of a universal quantum 
gate array on a quantum computer can only be realized in a probabilistic fashion, not in 
a deterministic fashion. By extending the one-to-one processor to one-input-to multiple- 
output processor, Yu et al proposed an approximate and probabilistic programmable multi- 
output quantum processor 10]. Schuch et al proposed a general programmable network to 
implement arbitrary two-qubit phase-shift operation with a network of one-qubit and two- 



qubit operations 



ll| . It embodies the following transformation \x) ^ e *^Ha;)w 
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over all possible basis states. Other interesting results were also obtained 

In quantum computation, a special controlled phase gate is of particular importance: it 
only changes the phase for a particular basis state |s), and leaves the other basis states 
unchanged. For instance, in the Grover quantum search algorithm j^, there are two phase 
inversions, Ir which inverses the phase of the marked state and leaves alone the phases of 
other basis states, and Iq which acts likewise. Starting from this controlled phase gate, other 
controlled gate operation, such as the multi-qubit controlled-not gate can be realized, and 
multi-qubit controlled gate is the major operation in many practical application such as 



initializing a quantum computer from |00 • • -0) into an arbitrary superposed state [16[, and 
in designing quantum cloning machines |l7|. This multi-qubit gate can be decomposed into 
a set of one-qubit gate and two-qub it g ates as given in Ref.J7[. In fact, this decomposition 
has been used in, for instance, Refs.|l8l. l21| to implement multi-qubit controlled-NOT gate. 
In the same spirit, the multi-qubit controlled phase gate can be implemented systematically 
using the method proposed in Ref. |ll|. 

In this paper, we propose a different realization scheme for the n -qubit controlled phase 
gate using modular structure and recursion. The network retains the same form when the 
system extends to more qubits. It contains an n-qubit controlled phase rotation about the 
z-axis plus an n — 1- qubit controlled phase gate. Again the n — 1- qubit controlled phase 
gate can then be further realized similarly. In one hand, the controlled phase gate is used in 
many computing algorithms, and it merits special attention. In the other hand, the modular 
realization of this controlled phase gate may provide convenience for some physical systems 
such as the nuclear magnetic resonance (NMR). Das et al proposed the scheme to implement 
the controlled phase-shift gate in the NMR system using line-selective pulses, where a con- 
trolled phase rotation was realized through three line selective pulses jl9|]. A line-selective 
pulse is the high- selective radio-frequency pulse designed to perturbed transverse magneti- 
zation one line at a time [20|. Experimentally, line-selective pulses can simplify a network 
circuit greatly, and makes it easy to extend the circuit to more qubits. The operational time 
of a line-selective pulse is not influenced by the number of the qubits but the spin system 
with well resolved couplings J2l|. For some certain systems, the line-selective pulses can 
use short time and this is an attractive feature considering the limited decoherence time in 
candidate quantum computer realization systems. Das et al demonstrated their scheme in 
a three qubit NMR system by simulations and in a two qubit system experimentally. Their 
method can be applicable for both the weakly-coupled spin systems and the strongly coupled 
systems. 

Our work concentrates on the weakly-coupled spin systems, which have been widely used 
in quantum computation J22]. We have experimentally demonstrated the feasibility of the 
proposed network in a three qubit NMR quantum system, where the multi-qubit controlled 
phase rotation is realized by the line-selective radio-frequency pulses and the spin selective 
pulses on the target qubit. Compared with the method used in Das et al's work, our method 
is experimentally simple. In our experiment, the line-selective pulse that pin-points directly 



the right frequency for the controlled configuration has been successfully demonstrated. The 
line-selective pulse has the advantage of the reduced operating time and ease in algorithm 
design. 

The paper is organized as follows. In sectional we describe the modularized realization of 
the multi-qubit controlled phase gate. In section ITTTl we report the details of the experimental 
realization in a three-qubit NMR quantum system. In section IIVI we discuss some issues in 
the realization and a briefly summarize the result. 

II. MODULARIZED QUANTUM NETWORK FOR CONTROLLED PHASE 
GATE 

The controlled phase- shift gate is one of most important operation in quantum compu- 
tation. On one hand, with a two-bit phase gate and one-qubit rotation gate, one can build 
any arbitrary unitary operation in quantum computation. On the other hand, multi-qubit 
controlled phase gate is the essential operation in many quantum algorithms, such as the 



Grover quantum search algorithm, quantum Fourier transform and so on[2, |23|,|2J,|25|. The 
controlled phase- shift operation for ra-qubit system is denoted by Iu\, where |s„) is a basis 
state of an ra-qubit system. If,^ transforms |s„) to e^'^\sn), and leaves the other basis states 
unaltered. Without loss of generality, we assume that qubit n is the target bit, and the 
other n — 1 qubits are the controlling bits. Let \sn) = |1 ■ ■ ■ 1), denoting that all qubits lie 
in state |1). Ify can be expressed as 

/ 



TV 



u ■ 



where 



e'^f 

and 12^-2 denotes the (2" — 2) x (2" —2) unit matrix. U is an operation applied to the target 
qubit. U cannot be constructed from the Pauli operators ax, cfy, cr^ and I2 directly, because 
it is not in SU{2) group[7j. Practically, U can be realized in the following manner, noting 
that 



e^'^/2 
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. . <l>(^/2)i?,(-v.), (3) 
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where Rz{—ip) = e" 



-iiph 



Iz is z-component of the angular momentum of the spin, and 



by a rotation about 
. For a single qubit, 



-^^lO) = ||0), /2II) = — 1|1), by setting Ti = \. Rz{—ip) can be realizec 

the 2;-axis, and in NMR it is realized by radio-frequency (rf) pulses (26 

the difference between Rz{—(f) and U is negligible because they differ only by a global 

phase. But for the controlled phase gate, this difference cannot be ignored. For instance, in 

a two-qubit system, the controlled- rotation R^{<f) is 
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and it is obvious the overall phase between U and Rz is no longer trivial, and ^{ip/2) must 
be considered in designing the network for If, y 

The network to implement If, ■. is shown by Fig. ^ where the "global" controlled-$((y9/2) 
is explicitly drawn. By a simple observation, it is not difficult to show that the controlled- 
$((/)/2) is equivalent to: 
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where I\s^_^\ denotes the phase- shift operation for the n — 1- qubit system composed of the 
control qubits in the original n-qubit system. |s„_i) is the basis state of the first n — 1 control 
qubits where the corresponding basis state of the whole n qubit is |s„) = |s„_i)|l)„. The 
ra — 1-th qubit is the target qubit, the other n — 2 qubits are the control qubits in If _ y The 
effect of this shortened controlled phase gate can be seen as follows: the controlled-phase 
rotation R'^{f) changes 



|Sn-lO) 
k„_ll) 



e-*'^/2|s„_i0). 



Jip/2\ 
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and it leaves other basis states unaltered. A follow-up I'fJ^_^\ operation changes |s„_i) 
into e''f''^\sn-i), and this takes e-''f''^\sn-i^) into |s„_iO), and e*'^/2|s„_il) into e*'^|s„_il). 
Hence the total effect is to change |s„) into e*'^|s„). Fig. [T] can be further represented 
as Fig. El where I\J^_^\ is applied to the control qubits. We call the I\J^_^\ part as the 
compensatory network module. The controlled-phase rotation can be realized by two multi- 
qubit controUed-NOT gates (A„_i(NOT)) and two single-bit phase rotations |7|. A„_i(NOT) 



can be directly realized in some physical realizations such as the NMR system, or be further 
reduced to a set of single qubit rotations and two-qubit controlled NOT gates as described 

f.Q. 



in Ref . 
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The nice feature in FigO is that J,'^ > can be realized in the same way as If y it 
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can be realized by a compensatory module composed of the n — 2 qubits together with an 
n — 2 qubits controlled phase rotation, and the compensatory module unit can be similarly 
constructed shown as the right network in Fig. |21 Hence we have obtained a recursive 
way to construct the multi-qubit phase shift gate. In summary, If^ is constructed by the 
n-1 -qubit controlled phase rotation and the compensatory network module I'fJ^_-^\- Ifx is 
realized recursively. 



III. NMR REALIZATION OF THE MODULARIZED NETWORK 

In the experiment, we use a sample of Carbon-13 labelled trichloroethylene (TCE) dis- 
solved in d-chloroform. Data are taken at controlled temperature (22°) with a Bruker DRX 
500 MHz spectrometer. ^H is denoted as qubit 3, the ^^C directly connecting to ^H is de- 
noted as qubit 2, and the other ^'^C is denoted as qubit 1. The three qubits are denoted as 
CI, C2 and H3. By setting h = 1, the Hamitonian of the three-qubit system is 

H = -2tivJI - 2%V2ll - 27rz/34^ + 2tx J^^llll + '2ixJ2zllll + 27r Jl3/i4^ (7) 

where z/i, V2) ^s are the resonance frequencies for CI, C2 and H3, respectively, and ui = 
z/2 + 904.4Hz. The coupling constants are measured to be Ju = 103.1 Hz, J23 = 203.8 Hz, 
and Ji3 = 9.16 Hz respectively. 

We realize Jm^^x in the two-qubit system composed of CI and C2, through decoupling H3. 
The spin-coupling evolution between CI and C2 is described by 

[r] = e-^2--^"^'^'. (8) 

CI and C2 construct a homonuclear system, which can demonstrate the implementation of 
If-^^^ clearly, and represent the effect of the compensatory module. Fig. |21 is simplified to 
a two-qubit network, where CI is the control qubit, and C2 is the target qubit. The com- 
pensatory module is expressed as Ifl = ^^{(p/A)Rl{—(p/2), and is equivalent to R\{~ip/2). 
$^(79/4) only contributes an irrelevant overall phase factor before Ifn)- A„_i(A^OT) is the 



ordinary CNOT operation. The operation sequence CNOT-i?^((y9/2)-CNOT can be realized 



by 
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Hence K^-, is realized by the following pulse sequence 
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where (p is taken as the menus value during the experiment. 

The experiment starts with the pseudo-pure state prepared by spatial averaging. The 



following pulse sequence 
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transforms the system from equilibrium state to the pseudo-pure state 
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which is equivalent to state |00). The symbol 1/4 J12 denotes the evolution caused by the 
magnetic field for l/4Ji2 when the pulses are switched off. After the application of the 
Hadamard transform, |00) is transform to the superposition of states 



b) = (|00) + |01) + |10) + |ll))/2. 
It is obvious Ifii\\p) = (1 00) + |01) + |10) + e*'^|ll))/2, represented by density matrix 
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In Eq. flllll . the matrix elements p/(l,3) and p/(2,4) can be directly observed through the 
spectra of CI, and p/(l,2), p/(3,4) can be directly observed through the spectra of C2. 
The Hadamard transform is realized by [— f ]^'^ — [^Yx"^- Ifu) is realized by 
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Figs. El^a-d) show the results when —f/2 = 7r/4, 7r/2, 37r/4 and 7r, respectively. In each 
figure, the frequency centers of the doublets of CI and C2 are 124.14ppm and 116.91ppm, 
respectively. The right and left peaks in the doublet of CI correspond to p/(l, 3) and p/(2, 4) 
in Eq.(fTT|). and the right and left ones in the doublet of of C2 correspond to p/(l,2) and 



p/(3,4) in Eq.()lip respectively. The phases of the right peaks of CI and C2 hardly change 
with —^/2. They are chosen as the reference phases to observe the phases of the other 
peaks so that the phases of signals are meaningful |29|. The phase of the left peak (denoted 
by Ph2) changes via —^/l proportionably, which are shown in Fig. EJ 

To show the effect of the compensatory module, we add R\{ip/2) at the end of network 
shown in Fig. |21 to cancel the compensatory network module. The network without the 
compensatory network module implements the transformation 

/in n n \ 
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It is obvious that If-^-^^R^i^ip /2)\p) 
represented as 
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From Eq. ()13|1 . the signals of CI and C2 have phase difference, which can be observed in 
this homonuclear system. 
R\{({)/2) can be realized by 
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n theory, ip/2 = 4(z/i — z/2)27r(t + tp/2), where tp is the width of the n pulse selective for C2 
26[. However, the dependence of ip/2 on t needs be measured due to the errors caused by 
decoherence and imperfection in the pulses. Figs. Efa-b) show the results of I^-^xRl{ip/2)\p) 
with —1^/2 = 7r/2, vr. Compared with Figs. IH^b) and (d), consequently the phases in peaks 
CI and C2 are no longer the same. Theoretically this phase difference should be <f/2, and 
our experiment results agree with this theoretical expectation. 



IV. DISCUSSION AND SUMMARY 

We realize the network shown in Fig. IHl in a three qubit system by switching off the 
decoupling pulses for if 3. The network consists of two rotation operations for C2 and two 



8 



controUed-controUed-NOT gates (CCNOT or Toffoli gate). C2 is the target qubit, CI and 
H3 are the control qubits with control condition |0)i|l)3. The network implements the 



transform /loo'lf /|o(?> 



, which can be transformed to -^^iqh) 



by multiplying /|ooV|oii>- ^^ ^^^t, 
-^1001) -^1011) ^^ ^^^ controlled phase- shift operation for the system composed of CI and H3, 
and it is represented as 
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with the basis state order |0)i|0)3, |0)i|l)3, |l)i|0)3, |l)i|l)3. -^10)111)3 can realized by a similar 



way in Sec. IIIII /, 



Mooi) -^loii) transforms the superposition of states 



b)3 = (|000) - |001) - |010) + |011))/2 
to (|000) - e-*^/2|ooi) _ |010) + e*'^/2|oii))/2, expressed as the density matrix 
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where \p)-j, is obtained by applying [f]^'^ to |000). The elements (1,3) and (2,4) in the matrix 
can be observed in the spectrum of C2 directly. 

The pseudo-pure state |000) can be represented by 



Peff = (\l + Il){\l + Il){\l + Il) 



(17) 



Peff can be expressed as pe// =Pi+P2 + |-^8x8, where 



Pi = {\ll + \ll + llll + \)ll 



(18) 



P^ = \^^\ll^l\l% (19) 

and the last term |/8x8 can be ignored [30|. Quantum computation can start with pi and 
P2, respectively. The average of the two results is equivalent to the result obtained from 
the initial state Pe//- Because the results of applying [f ]^'^ to p\ have no observable signals 
in the carbon spectra, it is sufficient to use p2 as the initial state to observe the results of 
applying /|ooi\ -^loii) ^^ \v)'i- This fact can simplify the process of experiments greatly. 



\TL\'i _\r,^r,A^ _ [5:11,2 1 rJ3 



The initial state described as Eq. (fT^ is realized by [|]^ — \grad\z — [f J^'" — gj [vrj^ — 

sJi; - \Af^' - si; - \-< - sir; - {-Af^' - [-f]^'' - t^™^]- CCNOT gate can be 
approximately realized by a line selective vr pulse with frequency v^ + J12/2 — J23/2, which 
causes transition |001) ^^ |011). Figs. [3 show the results of /|ooi> -^|oii)l^)3- -'■^ ^^^ figure, 
the right peak, with frequency v^ + J12/2 + ^23/2, corresponds the element (1,3), and the left 
peak, with frequency Z/2 + J12/2 — J23/2, corresponds the element (2,4), which changes versus 
(/?/2. Figs. [7|(a-f) show the results of ^/2 = 0, 25.2°, 38.2°, 90°, 141.8°, 154.7°, respectively. 
The graph in Fig. |H1 shows the phase of the right peak (denoted by P/i3) versus '^jl. It can 
be fitted as a line with slope approximate to 1.95. In theory, slope is -2. The difference of 
sign results from the approximation of CCNOT gates. 

The modularized quantum network to implement the controlled phase- shift gate is pro- 
posed, and realized on a three qubit NMR quantum computer. The n- qubit network nests 
the n — 1- qubit network. This idea represents the essential thoughts of computer science. 
The quantum modules provide a systemic method to building practical quantum comput- 
ers. It is sure that the modularized quantum networks will play central role for large- scale 
quantum computer. 
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FIG. 2: The network in Fig. ^ which is represented by a multiple qubit-controlled phase rotation 



operation and a compensatory network module, denoted by /,' 
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which is equivalent to the 



multiple qubit-controlled <l>((/?/2) in Fig. ^ The compensatory network module can be further 
decomposed as the network circuit denoted by the dashed rectangle in the right network. 



13 



{ 



124 122 120 118 116 ppm 



jl 



124 122 120 118 116 ppm 



u 



J 



124 122 120 118 116 ppm 



u 



124 122 120 118 116 Ppm 



FIG. 3: ^^C spectra of trichloroethylene after If > is applied to the uniform superposition of basis 
states obtained through applying Hadmamard transform to [00). ^H has been decoupled. The 
phases of the right peaks of in the doublets hardly change. The phases of the left peaks, denoted 
by Ph2 (2 for the two qubit system), change with —if/ 2 proportionably. Figs, (a-d) are the spectra 

for —95/2 = 7r/4, 7r/2, 37r/4 and vr, respectively. 
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FIG. 4: The phase of the left peak in the doublet of CI or C2, which is denoted as Ph2, versus 
—i-pjl. The graph can be fitted as a line with slope equal to 1.96. The theoretical expectation is 2. 
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FIG. 5: ^^C spectra after /|'^-^>i?^(99/2) is applying to the uniform superposition of basis states, 
when —^jl = '7r/2, vr, shown by Figs, (a) and (b), respectively. 
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FIG. 6: The quantum network to implement Lqq(\ ^igiiv 
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FIG. 7: Spectra of C2 after /|oifif -^|on) is applied to (|000) - |001) - |010) + |011))/2. Only two 
NMR peaks appear in each spectrum. The phase of the right peak hardly changes. The phase of 
the left peak changes with 93/2, proportionably. Figs, (a-f) show the spectra for c/9/2 = 0, 25.2°, 
38.2°, 90°, 141.8°, 154.7°, respectively. 
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FIG. 8: The phase of the left peak in the spectrum shown in Fig. [3 which is denoted as Ph3, 
versus (p/2. The graph can be fitted as a hne with slope approximate to 1.95. In theory, slope is 
-2. The difference of sign results from the approximation of CCNOT gates. 
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